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Abstract 

A number of synchronous moons are thought to harbor water oceans beneath their 
outer ice shehs. A subsurface ocean frictionahy decouples the shell from the interior. 
This has led to proposals that a weak tidal or atmospheric torque might cause the 
shell to rotate differentially with respect to the synchronously rotating interior. 
Applications along these lines have been made to Europa and Titan. However, the 
shell is coupled to the ocean by an elastic torque. As a result of centrifugal and 
tidal forces, the ocean would assume an ellipsoidal shape with its long axis aligned 
toward the parent planet. Any displacement of the shell away from its equilibrium 
position would induce strains thereby increasing its elastic energy and giving rise to 
an elastic restoring torque. In the investigation reported on here, the elastic torque 
is compared with the tidal torque acting on Europa and the atmospheric torque 
acting on Titan. 

Regarding Europa, it is shown that the tidal torque is far too weak to produce 
stresses that could fracture the ice shell, thus refuting an idea that has been widely 
advocated. Instead, it is suggested that the cracks arise from time-dependent stresses 
due to non-hydrostatic gravity anomalies from tidally driven, episodic convection 
in the satellite's interior. 

Two years of Cassini RADAR observations of Titan's surface have been inter- 
preted as implying an angular displacement of ~ 0.24 degrees relative to syn- 
chronous rotation. Compatibility of the amplitude and phase of the observed non- 
synchronous rotation with estimates of the atmospheric torque requires that Titan's 
shell be decoupled from its interior. We find that the elastic torque balances the 
seasonal atmospheric torque at an angular displacement < 0.05 degrees, effectively 
coupling the shell to the interior. Moreover, if Titan's surface were spinning faster 
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than synchronous, the tidal torque tending to restore synchronous rotation would 
almost certainly be larger than the atmospheric torque. There must either be a 
problem with the interpretation of the radar observations, or with our basic under- 
standing of Titan's atmosphere and/or interior. 
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1 Introduction 



Synchronous, or near synchronous, spinning satellites with ice shells and sub- 
surface oceans are an intriguing class of solar system bodies. Part of their ap- 
peal is the possibility that the oceans could support life. Currently Jupiter's 
satellite Europa and Saturn's satellite Titan are the prime candidates for mem- 
bership in this class. 

The induced magnetic field of Europa as measured by the magnetometer on the 
Galileo spacecraft essentially proves the presence of a current- carrying, near- 
surface liquid layer, i.e. a salty subsurface ocean (IKivelson et al.l . 120001 ). All 
mechanisms proposed for creating the observed morphologies o f cracks on Eu 



ropa' s surface require the pre s ence of a nea r -surfa ce fluid layer ( IGreeley et al 



19981 : iGreenberg et al.l . Il999l : iHoppa et al.l . Il999l ). Most mechanisms also in- 



voke an immeasurably slow, super-synchronous spin of the ice shell due to 
tidal torques associated with Europa's orbital eccentricity. The eccentricity, 
currently 0. 01, is continually excited by resonant interactio n s with lo and 



Ganym ede ( jPeale et al.l . ll979l : l0jakangas and Stevensonl . ll986l ). IGreenberg and Weidenschilling 
( Il984j ) s uggest that tidal torq ues could drive a super-synchronous rotation of 
Europa. iGeissler et al.l (119981 ) argue that the presence of a subsurface ocean 
might allow the surface ice shell to rotate non-synchronously even if the core 
is locked in synchronous rotation. 

Models of Titan's thermal evolution are consistent with the presence of a 
subsurface ocean. Moreover, motions of Titan's surface features as tracked 
by the Cassini spacecraft' s radar have beeri interpreted as implying a slow , 
non- synchronous rotation fistiles et al.l . l2008ll2009[l2010h . lLorenz et al.l fl2008[ ) 
propose that the non-synchronous rotation is driven by a seasonally changing 
atmospheric torque acting on an ice shell decoupled from the interior by a 
subsurface ocean. 



The plan of our paper is as follows. As a result of centrifugal and tidal forces, 
the interiors of synchronous satellites have ellipsoidal figures whose longest 
axes point toward the parent planet. We derive these forces and resultant 
equilibrium shape in §2. In §3, we evaluate the increase in elastic and gravi- 
tational energy associated with the angular displacement of an elastic outer 
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shell away from its equilibrium position and derive the elastic torque. In §4, 
we outline a derivation of the tidal torque exerted on a synchronous shell due 
to the satellite's orbital eccentricity. We derive the angular displacements at 
which the elastic torque balances the tidal torque acting on Europa's shell and 
the atmospheric torque acting on Titan's shell in §5. The results of this exer- 
cise imply that the long cracks on Europa do not result from non-synchronous 
rotation of the ice shell driven by the tidal torque. They also cast doubt on 
claims of non-synchronous rotation for Titan. We devote §6 to an explanation 
for the formation of large cracks on Europa's surface due to geoid variations 
associated with episodic mantle convection driven by tidal heating. We end 
with a short discussion in §7. The main points of our investigation can be 
appreciated with order of magnitude estimates. Consequently more detailed 
considerations are relegated to a series of Appendices. 



2 Potential theory 



This section contains background material selected for its relev ance to our in- 



vestig ation. For details and derivations, the reader might consult iMurray and Dermott 



(1200(1 ). 



Consider a synchronously spinning satellite with an elastic outer shell, an 
underlying ocean, and a solid core. Our focus is on the elastic shell which has 
mean radius R and thickness d « R. For clarity, the figures and Appendix 
A. 8 refer to a simpler model of an elastic shell and homogeneous fluid interior. 
Suppose the satellite moves on a circular orbit. Then the sum of the rotational 
and external tidal potential it feels is given by 

Vapplied = f -^P2,0(COS 9) - P2,2(cOS 9) COS 20j , (1) 



where n is the orbital frequency, P2fi{,x) = (3x^ — 1 )/2 and P2,2{x) = 3(1 — x^). 
To evaluate the satellite's shape under \4ppiied, account must be taken of the 
potential, V^eif, due to the distorted satellite. Up to quadrupole order, for 
r > R, 




It is customary to define 

J2 = ~C2fi ■ (3) 
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Expressed in terms of the principal moments of inertia, 
2C-A-B , ^ B-A 



J2 



and 



2,2 



(4) 





to planet 

Fig. 1. Left: Satellite with subsurface ocean. The outer shell (unfilled region) is 
elastic and the interior (gray region) is a fluid of uniform density. The black dots 
are points of reference fixed to the shell. Right: Satellite shape is defined by the 
deviation, ^, from the mean radius, R. 

Next we assume that the total potential, Kotai = Vseu + Kppiicd is constant on 
the satellite's surface. This is appropriate if the deviatoric stress in the elastic 
shell vanishes. Expressing the shape of the surface hj R + C,{0,(j)) (Figure 1) 
and defining the parameter q = n^R^ /{GM), we arrive at 



R 



5q 



+ C'2,o) A,o(cosi 



+ C2,2KP2,2(cos0)cos20. (5) 



It is common practice to express the gravitational coefficients 6*2,0 and 6*2,2 in 
terms of the Love number k as 

^2,0 = -^, C2,2 = ^. (6) 
D 4 



The parameter k depends upon the radial variation of density and elasticity. 
For bodies in hydrostatic equilibrium, only the run of density matters. The 
fluid Love number is denoted by kf (see Appendix A.l). Assuming k = kf, we 
obtain 
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3 Energies of Orientation and the Elastic Torque 



Consider a shell whose initial state is one of synchronous rotation and van- 
ishing deviatoric stressQ Suppose a spin-aligned torque applied to the shell 
rotates it by a small angle A << 1. Both elastic and gravitational energies 
increase. Two limiting cases are illustrated in Figure 2. In the first, the shell 
slides over the ocean which maintains its equilibrium shape. In the second, the 
entire satellite rotates with its shape remaining fixedH] In the simple example 
of a satellite with uniform density, only the elastic energy increases in the first 
case, whereas in the second only the gravitational energy does. 



31 « 1 



31 » 1 





Fig. 2. Two limiting cases for a satellite rotated away from its equilibrium position. 
Tlie line indicates tlie long axis of the ellipsoidal figure. If 7?. << 1, the shell deforms 
and stores elastic energy and the gravitational energy does not change. If 7^ >> 1, 
the shell behaves rigidly and the satellite's gravitational energy increases. Europa 
and Titan are both in the former category, TZ « 1 (see Table 1). 

In the first case, the meridians of the shell are stretched and compressed due 
to the rotation giving rise to a strain ~ Xq. The elastic energy of the shell after 
this rotation is i^eias ~ ^X'^q^dB? where is the shell's rigidity. In the second 
case, the equilibrium figure with a quadrupole [B — A) ~ qMR^ is rotated 
against the gravitational torque of the planet, increasing the gravitational 
energy Eg^av ~ qMR^n'^X^ ~ {GMyR)q^X\ 

The ratio of the energies in first and the second cases is 7^ = E^xas/E^ 



grav 



^ The assumption of vanishing deviatoric stress might be challenged. This point is 
discussed in §7. 

^ In order for the shell not to touch the solid core, the ocean must be ~ qR ~ 
few X 10^ m or greater. Titan's and Europa's oceans are estimated to be ~ 100 — 200 
km deep (jSohl et al.l . l2003l : lAnderson etaP . [l998,)- 
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fidR^ / {G M'^) . Next we evaluate the numerical coefficients of these energies 
without assuming uniform density. 



Expressions for the stress and strain tensors, cXjj and provided in 

Appendix A. 2. The elastic energy density generated by rotating the shell 
an angle A with respect to the equilibrium figure of the ocean is 

S = aijUij/2. (8) 

The total elastic energy 



Eelas — dR? j dilS 



487T fl + v 



;i + kffq'^idR^X^ 



(9) 



with V Poisson's ratio. 

The increase in gravitational energy due to rotating the entire satellite with 
its shape fixed so that its long axis is displaced by an angle A from the 
direction toward the planet is 



grav 



3 GM2 

-kf q X 

2 ^ R ^ 



(10) 



The ratio 



n 



E, 



clas 



E^ 



grav 



327r /I + z/ 
5 + u 



[l + kfyfidR^ 
Ff GAP 



In a real rotation the total increase of energy would be minimal. Thus if 
TZ << 1, the shape of the shell would conform to that of the ocean's equilib- 
rium figure, whereas if 7^ >> 1, the shell would maintain its shape and the 
ocean's figure would ro t ate with it. For Europa and Titan, 7^ << 1 (Table 1) . 
Karatekin et all fl2008h : lvan Hoolst et all fl2008[ ): and Ivan Hoolst et al.l fl2009h 
implicitly assume that the shell has infinite rigidity, or 7^ >> 1. Henceforth 
we take 71 « 1. 



The torque associated with the storage of elastic energy follows as 



dE, 



elas 



elas ■ 



967r /I 



5 + 1/ 



;i + kfYq'^idR'X 



:i2) 



A massless membrane provides the simplest model for a shell. A rotated and 
deformed membrane exerts a pressure on the ocean's surface due to the hoop 
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stress, 



Pi = {o-e,e + cr^,^) 



R 



-12 



.2 e sin 20 . 



(13) 



The elastic torque arises from the back action of the ocean on the shell as 
indicated by Equation A13. The primary effect on the ocean is to rotate the 
long axis of its ellipsoidal figure by an angle 5 ~ T^-AH] In the simple case of 
a massless membrane, the planet's gravitational torque balances the pressure 
torque the shell exerts on the ocean. We evaluate 6 for a more realistic model 
in Appendix A. 7. 

The surface density perturbation that results from the deformation of the shell 
gives rise to an additional pressure torque —{pgR/2fi)Teias acting on the shell 
as shown in Equation A170 This factor multiplying Teias is of order unity 
for Europa and Titan. The minus sign follows because hoop stress is maximal 
where the shell is stretched thinnest. It is noteworthy that the pressure torque 
on the shell due to the surface density perturbation is exactly canceled by 
the sum of the gravitational torques from the satellite and the planet on the 
distorted shell; we refer the motivated reader to Appendix A. 6 for a complete 
derivation. 



4 Tidal Torque on the Shell 

Tidal torques are often invoked to explain the origin of the putative non- 
synchronous rotation of Europa's ice shell. A tidal torque acts on the syn- 
chronously rotating shell of a satellite that moves on an orbit of low eccen- 
tricity. The torque arises from energy dissipated by the strain produced in the 
shell from diurnal tides. A tidal wave can travel a distance R across an ocean 
of depth if in a time R/yfgH. This time is much shorter than the diurnal 
time of Europa and Titan for an ocean depth of 100 km. The tidal response 
of the ocean is essentially instantaneous. The free libration frequency of the 
shell (with restoring torque due to elasticity) is an order of magnitude less 
than the diurnal frequency, n. Therefore to a good approximation, the shell 
remains stationary in the rotating frame and is deformed by the ocean tide be- 
low. Time dependent components of the diurnal tide can be sorted into terms 
with inertial pattern speeds Vtp = (n/2, 3n/2). Stresses and strains in the shell 
due to diurnal tides are given in Appendix A. 2. Corresponding elastic energies 
stored in the shell follow directly. 



^ The change in ellipticity is of second order in 5 and is ignored. 
The subscript i denotes ice. 
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En/2 = ^vr {h,eq)'fidR' , (14) 

Esn/2 = ^vr (/i,eg) W , (15) 

where ht is the tidal Love number. We expect ht < hf since elastic stresses 
from diurnal tides do not have time to relax. A mple expression for ht is given 
in Equation A3. The torque exerted on the shell by each component tide is 
related to its rate of energy dissipation hj T = E/ (Qp — n), where E = 2nE/Q 
follows from the definition of Q as 27r times the peak stored energy divided by 
the energy dissipated over a cycle. Each component tide has frequency n in 
the synchronous frame and the stored elastic energies are half the peak energy 
densities multiplied by the volume of the shell. Thus T„/2 = — 4£'n/2/<5, and 
^3n/2 = 4:E3n/2/Q, giving the total tidal torque 

^ 576 /I + z/\ 

Te-tide = — TT —— -^{qef^dR" . 16 
5 \b + uj Q 



The tidal torque tends to increase the spin rate because the tides are strongest 
near periapse where the orbital angular velocity exceeds the mean motion. 



5 Applications: Europa and Titan 



We now apply the results derived in earlier sections to Europa and Titan. 
We adopt parameter values /i = 4 x 10^° dyne cm~^, z/ = 1/3 and Y = 10^ 
dyne cm~^ for ice. The yield stress of ice, Y, is the most crucial and also the 
least well-constrained. We assume that the interiors of Europa and Titan have 
relaxed to hydrostatic equilibrium. Models predict kf ^ 1 for both bodies, a 
value that i s consistent with gravi ty measurements made during Galileo flybys 



of Europa (lA nderson et al.l . 119981 ) . 



We make use of four derived quantities. The first two, q and TZ, were defined 
previously. The third, the angle Aext at which the elastic torque. Equation 12, 
balances an externally applied torque. Text, is 

_ 5 f 5 + u \ Text . . 

9671 + u J {1 + kfyq^fidR^ ' ^ ' 



In the special case Text = Ttidai, 

2 



_ 6e^ / ht 
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which is much smaller than the angle of geometrical libration. In this case the 
thickness of the shell drops out. The fourth quantity is the forced libration 
amplitude 



7o 



9 2 ' 



(19) 



with eccentricity e and (elastic) free libration frequency ul = Teias/C*sheiilZ] 
The moment of inertia of the shell Csheii = {8tt /3)pR^d with p the density of 
ice. 



Parameters for Europa and Titan are given in Table 1. 
Table 1 

Parameters for Europa and Titan. 





Europa 


Titan 


M 


4.8 X 10^5 g 


1.4 X 10^6 g 


R 


1.6 X 10^ cm 


2.6 X 10^ cm 


e 


0.01 


0.03 


n 


2 X 10"^ 


4.6 X 10"^ s"^ 


Q 


5.4 X 10"^ 


3.9 X 10"^ 


n 


2.1 X 10"2((i/10km) 


0.1(d/100km) 


At 


2.3 X 10^^/Q rad 


2 X l{)-^/Q rad 


Aa 


N/A 


7 X 10"''((i/100km) rad 


7o 


52" 


5" 



5. 1 Europa 



IZ « 1 for Europa so to a good approximation its long axis points toward 
Jupiter. The elastic torque balances the tidal torque at an angle At ~ 2.3 x 
10~^/Q radians. We have assumed an ht defined in Equation A3. Rotation 
of the shell by this angle generates stresses of order qpXt ~ 5 x lO'^/Q dyne 
cm~^{£] These stresses are clearly inadequate to crack the shell even for Q = 1. 
The thickness of the elastic shell does not enter into this calculation. Since 
for Europa n > Uq, the shell does not follow the tide and the shell's forced 
libration amplitude 7o ~ 52". 



^ This amplitude and frequency is relevant for satellites with TZ « 1. 
^ The maximum stress is a factor of 4.5 times larger. 
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5. 2 Titan 



Titan possesses a thick atmosphere whose wind patterns exhibit seasonal varia- 
tions resulting in the exchange of angular momentum with the surface. Model 
estimates yield AL = 1.3 x 10^^ g cm^ 



^ for the total seasona l angular 
2009h . The ex- 



2005; Mitchell 



momentum exchange (iTokano and Neubauer 
change occurs at twice Saturn's orbital frequency, w, giving rise to a maximum 
torque 



T'atmosphere = 2wAL ~ 2 X 10 dyne cm . 



(20) 



This torque is capable of forcing a 100 km thick shell to rotate by = 7x 10""^ 
radians relative to the equilibrium shape of the ocean's surface. This angle is 
small in comparison to the angular disp l aceme nts of 0.12 deg/year reported 
by the Cassini radar team (jStiles et al.l . |2010| ). Moreover, since the shell is 
effectively coupled to the interior, the angular displacement would be in phase 
with the torque instead of half a cycle out of phase as it would be if the shell 
could rotate independe ntly of the interior. This latter point was raised by 
van Hoolst et al.l ( 120091 ) , but they came to this correct conclusion by incorrect 
arguments for reasons presented in §3. Since for Titan n » Uo, the shell does 
not follow the tide and the shell's forced libration amplitude 70 ~ 5". 



There is an additional problem with the notion that Titan's surface is spinning 
faster than synchronous. If it were, tides raised by Saturn on Titan would 
contribute an additional torque ~ 3 x lO^^/Q dyne cm tending to restore 
synchronous spin. Unless Q > 10^, a remarkably large value for a solid body, 
the tidal torque would overwhelm the atmospheric torque. 



6 Mantle Convection, Gravity Anomalies, &z Cracks in Europa's 
Ice Shell 



Having established that non-synchronous rotation did not cause the long 
cracks observed in Europa's ice shell, our challenge is to come up with a more 
promising hypothesis. Here we propose that like lo, Europa undergoes episodic 
tidal heating resulting in time-dependent, non- hydrostatic stresses which per- 
turb its external gravitational potential. The shape of the ocean surface and 
ice shell would respond to changes in the gravitational potential with cracks 
in the shell being a consequence of these changes of shape. 

Two pieces of evidence point to episodic heating of lo. Estimates of the current 
heat flow from its interior, P ^ 10^^ watt, exceed the firm upper bound of 
3.5 X 10^^ watt on the time-averaged tidal heating rate of its interior (see §A.9). 
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A unique feature of episodic tidal heating is t hat it can cau s e a t emporary 
shrinking of lo's orbit, just what is found by iLainey et al.l (120091 ) in their 
recent analysis of the measured orbital motions of lo, Europa, and Ganymede. 

Europa presents a more problematic case for episodic heating than lo does. 
Observations do not offer evidence for either volcanism or an enhanced heat 
flow. However, the paucity of impact craters on Europa's surface suggests 
that it was resurfaced as recently as 10^ years ago, perhaps as the result 
of an episode of enhanced tidal heating. Europa has a larger orbit and a 
smaller mass and radius than lo. Each of these tends to reduce the rate at 
which tides dissipate energy in its interior relative to that of lo. Even so, 
based on a Q = 100 and a tidal rigidity scaled to that of Earth's Moon, 
Europa's eccentricity damping timescale is only 3 x 10^ years. Moreover, on 
average, Europa's orbital eccentricity is likely to be much larger than Jo's as 
it is pumped by Europa's two-body mean motion reson ances with th e more 
massive satellites lo and Ganymede. The findings from iLainey et al.l (120091 ) 
suggest that at times lo must be deeper in resonance with Europa than it 
currently is. At such times Europa's orbital eccentricity would be larger than 
its present value. 



T, 



AT=T2-Ti 



Tidal Heating 



▼ 

V 



t 

V 



R 



Fig. 3. Mantle convection. A conducting lid of depth 6 caps a convecting layer 
driven by negative buoyancy of the cooler lid. The total depth of the mantle is R. 
The temperature drops from T2 to Ti < T2 across the lid. Density anomalies in the 
mantle from the subducting slabs vary on the convective time R/v. 

Our proposal, although unavoidably speculative, is guided by the c ase for 
episodic tidal heating of lo made by lOjakangas and StevensonI (119861 ). They 
investigate the response of lo to tidal heating based on detailed models for 
its rheology. Their calculations include, the angular momentum transfer from 
Jupiter that drives the expansion of lo's orbit, the excitation of lo's orbital 
eccentricity arising from its mean motion resonance with Europa, and the 
damping of the eccentricity due to tides raised in the satellite by Jupiter. 
In plausible scenarios, the satellite undergoes episodes of enhanced heating 
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lasting of order 10^ years separated by quiescent intervals with durations of 
about 10^ years. Good matches are obtained for both the current values of 
lo's heat flow and orbital eccentricity. Now we investigate the consequences of 
assuming that Europa undergoes episodes of en h anced tidal heating similar 



to those proposed by lOjakangas and StevensonI (119861 ) for lo. Suppose just 



10^'^ watt were dissipated and stored inside Europa during an interval of 10^ 
yearsQ If spread uniformly throughout the satellite, this heat would raise 
its internal ternperature by AT ~ 50 K and cause its radius to expand by 



AR/R ~ 10"'^ 1^ The tensile stress associated with the expansion would be 
of order 4 x lO^dynecm"^, well above what what is required to crack ice. 
But tidal stresses are not spatially uniform so heterogeneous heating is a 
more realistic option. It might generate convective plumes similar to the one 
illustrated in Figure 3. These would generate non-hydrostatic perturbations 
of the gravitational potential over a wide range of scales. The satellite's shell 
would adjust to an equilibrium shape in which the elastic and gravitational 
forces were in balance. Next we provide an estimate of the scale dependent 
elastic stresses. 

Consider an elastic shell whose relaxed state conforms to the hydrostatic shape 
of the synchronously rotating satellite. The vertical displacement of the shell 
relative to the relaxed state is described by where i denotes angular degree. 
The displacement gives rise to strains that scale as 



1 + ^U^. (21) 




The first term in the parentheses arises from stretching and the second from 
bending. The latter is appropriate for a flat plate so its dependence on i might 
seem surprising. However R/i is just the arc length associated with an angular 
separation 1/^ on a sphere of radius R. The elastic energy is 

Ei,ei.s ~ fiR^dej ^id{l + ed/Rfij. (22) 



The gravitational energy associated with the departure of the ocean's surface 
from an equipotential is given by 

^£,grav~^pi?'(0-6)' , (23) 



where is the height of the geoid. The shell deforms to minimize the sum of 



^ About 10% of the current heat flow from lo. 

^ Europa's mantle accounts for ^90% of its radius ( Anderson et al. . 19981 ). 
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the elastic and gravitational energies. Thus 

^ Y+TV^l + Pd/Rf ' ^^^^ 



and 



(1 + fd/R)Q/R 



Here 7?. ^ 1 is defined in Equation 11. 



If the equipotential were to obey Kaula's rule (iKaulal . Il966l ). which approxi- 
mately characterizes the scale dependence of the surface gravitational poten- 
tials of Earth, Venus, Mars, and the Moon, the rms potential difference would 
vary in direct proportion to the inverse angular separation, A0 oc £-1.0 The 
height of the geoid would scale similarly, Q ~ ^'P/q ~ the largest 

scales ^£ ~ Q and the shape of the shell would conform to that of the surface 
equipotential giving rise to strains that decrease as ~ Cq^-^ with increas- 



ing ^12^ For ^ > (R/dy^"^, the elastic energy from bending the shell exceeds 
that from stretching. In this regime the strain would increase linearly with £, 
€£ ~ eo{d/R)i, up to a scale iu given by 



For i > Imi the vertical displacement of the ice shell would be significantly 
smaller than the height of the geoid. 



1 R^ 



7^£4 




(27) 



Thus the strain would be given by 

J_f!MV<i~i^^^ (2f 
' \ t R RdP' ■ 



The parameters listed in Table 1 for Europa yield R^/'^/d^/'^ ~ 130 km for 
the length scale at which the strain is minimal, and R£m ~ 50 km for the 



^ At least for Earth and Venus, the low order components of the non-hydrostatic 
part of the gravitational potential are thought to arise from mantle convection. 
Here eo is the strain at small I. 
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length scale at which the surface no longer conforms to the geoid. Global 
cracks are the primary focus of our paper. However, we note that the scale of 
Europa's cycloidal cracks is similar to that at which strains due to stretching 
and bending have comparable magnitudes. Moreover, an interplay between 
strains due to bending and stretching accounts for t he formation of cycloidal 
crack s when blunt objects tear thin, brittle sheets (iGhatak and Mahadevaru . 



20031 ). Perhaps Europa's cycloids form from the blunt tearing action of ~100 



km scale, time-dependent, geoid anomalies. 



7 Discussion 



We conclude with a few parting thoughts. 

Our calculation of the elastic torque follows from the assumption that the 
shell's deviatoric stress vanishes in its tidally deformed state. But there is an- 
other possibility!"! The shell might have formed rotating super-synchronously 
with a deviatoric stress that would vanish if it were to assume an axisymmetric 
shape. Then non-synchronous rotation could proceed without any variation of 
elastic energy and thus unimpeded by an elastic torque. Although this is a 
logical possibility, the weakness of the drivers of non-synchronous rotation for 
Europa and Titan render it implausible. The tiny values of At calculated in §5 
are a measure of just how axisymmetric the unstressed shells would have to 
be in order to sustain non-synchronous rotation. 

We have demonstrated that Europa's long cracks are not opened by stresses 
due to non-synchronous rotation driven by the torque arising from the di- 
urnal tide. Although this removes much of the evidence for non-synchronous 
rotation, it does not preclude it. Provided the small stresses estimated in 
§5 could relax on timescale treiax, nonsynchronous rotation might proceed on 
timescale treiax/At. Our use of might instead of would is deliberate. If the shell 
were to spin faster than synchronous, it would be subject to a component of 
tidal torque arising from the semi-diurnal tide which would tend to restore 
synchronous spin. In the limit of small orbital eccentricity, the semi-diurnal 
component of tidal torque would overwhelm the diurnal component unless the 
tidal Q were a steeply increasing function of period. 

Non-hydrostatic moment differences would drive the ice shell to reorient with 
respect to a system of axes pointing toward the planet, along the orbit, and 
normal to it. However, because the ellipsoidal shape of the ocean surface would 
maintain its orie ntation in these ax is, the elastic torque would oppose this type 



of polar wander ( IWillemanru . Il984j ) . Significant reorientation of the shell would 
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As pointed out to us by both Re'em Sari and Jay Melosh. 
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require that the ice either crack or creep. 

Provided Titan's shell is ~ 100km thick, the presence of an elastic torque 
presents a challenge to the interpretation of Titan's non-synchronous rotation. 
Observational evidence for the non-synchronous rotation of Titan cannot be 
lightly dismissed. Although we suspect an error in the data analysis, it is also 
possible that we are missing something important in our understanding of the 
properties of Titan's atmosphere and/or interior. 



A Appendix 



A.l Love numbers 



Under the action of an applied potential U, the surface of a self-gravitating 
body of radius R is lifted by ^ = hUg/g, where Ug is the external potential 
evaluated at the surface. The distortion of the body gives rise to an additional 
gravitational potential kUg- Love numbers h and k are derived below for simple 
cases relevant to this paper. 

The external potential acting on the satellite is taken to have angular depen- 
dence proportional to P2,o(cos7). The response of the satellite can be written 
as a surface density perturbation of the form E = EoP2,o(cos7). Its potential 
is given by 

Un = z . (Al) 



Example 1: An incompressible fluid body of density p. The surface is lifted by 
^ — H/p. Then from the definition of h, 

U 3 5 3 

^ = — -I- which implies h — - ^'k— - . (A2) 
g 5 2 2 



Example 2: A rigid sphere of density Pc and radius Rc overlain by an ocean of 
density p and surface radius R. Define p = [PcRl + p{R^ — -^c)]/-^^- Then 

^^^ + ^^ which implies /i=^^&k = ^^. (A3) 
g bp 5p — 3p bp — 3p 



Under loading by a layer of surface density S, the boundary condition at the 
top of the ocean becomes gT, — pU — const. Thus the distortion of the water's 
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surface satisfies 



-— + — + =eu, = — + ^ + A4 

gp 5g p p hp p 



which yields 



i^ = --. (A5) 
P 



A. 2 Stresses and Strains 



We follow the technique described in iMatsuyama and Nimmd (120081 ) to cal- 
culate components of the stress tensor associated with the reorientation of an 
elastic shell relative to the equilibrium shape of the ocean. 



A. 3 Stress From Rotation of Shell By Angle A 



After factoring out the coefficient 

(1 + kf)qp\ , (A6) 



we find 



aee = -T{2 + sm^e) sin 20 , 
= r(2 — 3 sin^ 9) sin 20 , 
(76)0 = 2T cos 9 cos 20 , 

c"</,e = o"e0- (A7) 



A. 4 Time Variable Stress From Diurnal Tide 



After factoring out the coefficient 

A = (l^) hteqp (A8) 

we obtain: for pattern speed n/2, 
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3 

cree = -^-^(2 + sin^ 9) cos(20 + nt) , 
3 

o'(t>(t> — ^^(2 ~ 3 sin^ 9) cos(20 + nt) , 
3 

ao4,^—-Acos9sm{2(f) + nt), (A9) 
and for pattern speed 3n/2, 
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(7^0 = -^A{2 + sin^ 61) cos(20 - nt) , 
21 

cr,^^ = — ^^(2 ~ 3 sin^ 61) cos(20 - nt) , 
21 

ag^^ — —A sin 9 cos{2(p - nt) . (AlO) 



A. 5 Strains from Stresses 

Components of the strain tensor follow from those of the stress according to 
2//(l + v) 

_ (^(jxj) — T^Crgg 

2/.(l + i.) ' 

U4,g = Ug^ . (All) 

A. 6 Pressure Torques 



The hoop stress from the deformed shell exerts a pressure perturbation, pi, 
on the ocean's surface; 

Pi = i.'Jee + 'J^xt')^ = -12 (^^) (l + A;/)gAi-^Asin2^sin20, (A12) 



with (Tgg, (70^ given in Equation A7 and values oikf derived in Appendix A.l. 
Acting on the shell, pi produces a torque Tp^ ; 

Tp,^-R^ Jdn^Pi, (A13) 
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where 



kf)q sm^esm2<j) (A14) 



and ^ is given in Equation 7. It is a simple exercise to demonstrate that 
Tp^ = Teias (Equation 12). 

We neglect the difference in the densities of water and ice, Pw = Pi = P- Pertur- 
bations in the surface density of the shell, S2 = —dpiuee + u^^) , from strains 
in Equation All are responsible for an additional pressure perturbation, p^] 

P2 = . (A15) 



Note that 

P2 f^-^\ 9PR 



(A16) 



so 



Since gpR/2p ~ 1 for Europa and Titan, the shell is more than a massless 
membrane and Tp^ cannot be neglected. 



The gravitational torque on the surface density perturbation S2 is 

T,, = -R' I dn^^E, , (A18) 

where 

^ = ^(1 + Mgsm ^sm20 = -,-, (A19) 



and Vtotai is the sum of Equation 1 and Equation 2. Tg^ exactly cancels Tp^ 
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A fraction 1/(1 + kj) of Tg^ is due to the planet and the remainder is due to the 
tidally distorted satellite. 
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A. 7 Distortion of the Surface 



The pressure perturbation, pi + p2 given in Equations A12 and A15, distorts 
the shape of the ocean's surface by a; 



a 
R 



htPi + P2 
pgR 



= 12 



1 — z/\ pgR 



l±y\ (1 + kf)q pd ^ ^.^2 
5 + i/y pgR^ 



Asin^^sin2(/) .(A20) 



Note that a is an order d/R smaller than ^, Equation 7. The long axis of the 
ellipsoidal figure rotates by an angle 



5_ 
A 



ht- 



1 



pgR 
1 + uJ 2ii 



1 + 1'\ pd 



5 + uJ pgR^ 



(A21) 



Two points are worth noting. 



• 5/ A is of order TZ which is small in applications we are considering. 

• The fractional change in ellipticity is of order S. 



A. 8 Tying It All Together 



We offer an overview of the torques acting on different parts of the satellite and 
show that they sum to zero. To maintain clarity, this is done in the context of a 
simple model. The satellite is taken to move on a circular orbits. It consists of 
an elastic ice shell and a fluid water interior. The deviatoric stresses vanish in 
its initial state. The density difference between unstressed ice and liquid water 
is neglected so A;/ = 3/2 and ht = 5/2. Imagine that a torque TappUcd = —Tdas 
(Equation 12) is applied to the clastic shell which rotates by a small angle, 
A, with respect to the equilibrium shape of the ocean that lies underneath it. 
The ocean's long axis rotates relative to the direction to the planet by the 
smaller angle S (Equation A21); 



5 
A 



20 



1 - 



2 pgR 



1 



pd 

5 + 1/7 pgR^ 



(A22) 



Stemming from its rotation by S, the planet exerts a torque on the ocean 



Ts = -3{B - A)n'^S 



2R 



1 - 



2 pgR 



5 \l + vj 2p 



elas ■ 



(A23) 
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It also exerts a torque on the perturbed surface density of the shell given by 



Summing these two torques confirms that the total planet torque equals Teias- 

The torques on the shell also must sum to Tdas- From the back reaction of the 
pressure it exerts on the ocean below, there is 

2pi — ^elas (A25) 



from the hoop stress (Equation A13), and 



from the perturbed weight of the shell (Equation A17). The gravitational 
torque includes 



from the planet acting on the shell's perturbed surface density (Equation A24), 
and 



3 (^^\ 9PRr 

5\l + iyJ 2fx 



Ii» = ?(^^)^^.,„ (A28) 



from the fluid interior acting on the shell's perturbed surface density. Summing 
the torques on the shell we have accounted for so far, we obtain 

-^shell ~l~ ~^ -^21 ~^ -^22 -^1 -^elas • (A29) 



A. 9 Bounding Tidal Dissipation Withm lo, Europa & Ganymede 



The angular momentum L and energy E oi a, satellite with mass M moving 
on a circular orbit with mean motion n around a planet of mass Mj may be 
written as 



n 



2/3 



n 



1/3 



and E = , (A30) 
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where /i = GMj. We express the tidal torque the planet exerts on a satellite 
by 

T = KM'^n'^ , (A31) 

with K is assumed to be independent of both M and n. 

Now consider three satellites, denoted 1, 2, 3 in order of increasing orbital 
semi-major axis. Conservation of angular momentum implies 

Li + L2 + L3 = Ti+T2 + T3. (A32) 

If tides raised within the satellites maintain their orbital eccentricities at low 
levels, the total (time-averaged) power dissipated within them is given by the 
difference between the rate at which work is done by tidal torques and that 
at which the orbital energy increases; 

P = n,T, + nsTa + n^n - ^1 - ^2 - ^3- (A33) 



Next we assume that the orbits are locked together by mean motion resonances 
such that rii = 2n2 = 4n3 and that we can neglect T2 and T3 with respect to 



Ti. on Combining equations (A33) and (A32) then yields 



P = KMln\ 1 



1 + 2-2/3(M2/Mi) + 2-^/^{M^/Mi 



[1 + 2V3(M2/Mi) + 22/3(M3/Mi)] 



(A34) 



Lissauer et al.l (Il984f ) derive a similar result for the heating of Enceladus from 
interactions with Janus and Saturn's rings. An upper bound on K is obtained 
by assuming that the resonance formed early and that subsequently the system 
expanded by a significant fraction of its current size during the age of the solar 
system, tgs- With these assumptions we find 



GM., 



ub 



l?,Miaitssn\ 



1 + 2^/^(M2/Mi) + 2^/3 (Ms /Mi; 



(A35) 



Consequently, 

2i/3G'Mj 



ub 



M2 + 



22/3 



(A36) 



With values appropriate to the Galilean satellites. Pub = 3.3 x 10 watt. 
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Both assumptions are valid for the Gahlean satelhtes. 
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